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Abstract 



In this paper, we study gamma ray spectra for various scenar- 
ios of dark matter annihilation and decay. We focus on processes 
which generate only high-energy photons or leptons and photons, but 
no proton-antiproton pairs, to be compatible with PAMELA'S data. 
We investigate photons produced directly from two-body decay chains 
and photons produced together with charged particles. For the former 
case we also include the process DM(+DM) — > Ncf) 2Nj which can 
arise from specific strongly-coupled dark matter scenarios. For the 
latter case, photons are either generated by final state radiation from 
high-energy leptons or are directly generated from contact interactions 
represented by higher-order (non-renormalizable) operators obtained 
after integrating out heavy modes. We compare their overall annihila- 
tion cross-sections/decay rates taking into account chiral suppression 
(in the s-wave approximation), dimension of operators and dark mat- 
ter particle properties. A rough estimate shows that, for a dark matter 
particle with a mass of O(lTeV), the hard photon spectra in direct 
electron-positron-photon final states arising from either scalar boson 
dark matter annihilation/decay or Majorana fermion dark matter an- 
nihilation are dominated by higher-order operators if the scale of the 
leading operator is lower than 0(1000 TeV). Otherwise, all the pho- 
ton spectra arising in this way are dominated by final state radiation. 
Among the spectra studied, the higher-order operators spectrum is the 
hardest while the final state radiation spectrum with an intermediate 
decay is the softest. 
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1 Introduction 



Recently several experiments (ATIC [2], H.E.S.S. [3lll], PAMELA [5], FERMI 
[6]) detected an anomaly in the cosmic ray electron spectrum, measuring an 
excess in the high-energy positron flux compare to usual diffusion models. 
Such an excess was not found for antiprotons [7]. One was thus led to conjec- 
ture a new source of primary electrons and positrons. The most likely sources 
advanced to explain this anomaly are nearby pulsars or dark matter (DM) 
annihilation/decay [El [9]. In any case, the issue related to the identity of 
the new source (pulsars versus DM) will be clarifled by the photon spectrum 
which will be measured by FERMI and announced later this summer. In- 
deed, there is an irreducible background of gamma ray photons coming from 
flnal state radiation (FSR) from electrons and/or positrons. Nearby pulsars 
would lead to a local photon spectrum while DM would instead generate a 
diffuse photon spectrum. 

The existence of DM is well established, although its paticle identity is 
still unknown. Assuming that DM annihilation/decay is the new source of 
primary electrons and positrons, several scenarios are possible (see [10-33] 
for related work). Apart from the FSR of photons, DM might also annihi- 
late/decay directly to photons. Different photon spectra are expected from 
each annihilation/decay mode and a knowledge of the achievable spectra will 
help physicists understand the properties of DM at the particle level. 

In Section O we discuss direct production of photons from DM annihi- 
lation/decay through subsequent two-body decay chain. In Section [3l we 
study the irreducible photon background from charged particles, consider- 
ing only the dominant process (either FSR of photons or photon production 
from higher-order operators). Next, in Section H] we consider direct photon 
production from taus. Finally, Section O contains a summary of the results 
and a comparison of the different photon spectra and total fluxes obtained 
from the density of states (DOS). Various detailed computations are left for 
the Appendix. 

Since we want to address both annihilation and decay, we compute pho- 
ton DOS throughout the paper, and it only differs with the spectrum by a 
constant factor as total number of photons. One can write DOS as 

1 dN^ 1 Mav) 1 dN^ 1 dT 
L = ; — L or = 

N^dE^ {av) dE^ dE^ T dE^' 

not photon spectra The photon DOS are thus normalized to 1, i.e. 

J dEry^^^ = 1, except for FSR where the DOS is normalized with respect 
to the annihilation cross-section/decay to zeroth-order in the flne structure 
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constant (see Section [3]). Thus, the total number of photons A'^^, given by 
= J dE^^^, can be found from the photon multiphcity. For example, 
a scenario, where DM annihilates/decays to two scalar bosons which subse- 
quently decay to two photons each, would have a total number of photons 
given by A^^ = 4. 

Since we focus on the photon DOS, all the results are applicable to DM 
annihilation and decay, and both are assessed simultaneously by substituting 
the parameter M by 2mDM for DM annihilation or for DM decay. Since 
backgrounds fall roughly like E'^, we plot the photon DOS in function of 
the dimensionless photon energy 2E^/M with an extra E'^ factor. This is 
consistent with the standard representation used in experiments. The M/2 
factor appearing in the plots is to make the quantities dimensionless, where 
M/2 is the maximal energy a photon can get during a process. 

Throughout this paper we assume that DM annihilates/decays only to 
leptons or photons as indicated by the experiments. Moreover, we assume 
that DM annihilation always occurs in the s-wave approximation. Under this 
assumption, we can get rid of the operators that can only contribute in p 
wave or higher order waves when we do operator analysis in later sections. 
Furthermore, DM annihilation to leptons or photons is allowed whatever the 
particle identity of DM (scalar boson, fermion or gauge boson). However, 
with the assumptions that standard model particles are not charged under 
hidden symmetries and that individual lepton numbers are conserved, DM 
decay to leptons or photons is allowed only for scalar boson and abelian gauge 
boson DM. 

Issues related to the overall annihilation cross-sections/decay rates (Som- 
merfeld enhancement, non-thermal DM production, etc) and the irreducible 
astrophysical photon background (inverse Compton scattering from starlight 
and CMB, synchrotron radiation from galactic magnetic fields) will not be 
investigated. 

2 Direct production of photons through sub- 
sequent two-body decay chain 

In this section, we analyze direct production of photons through subsequent 
two-body decay chain produced by DM annihilation/decay. We assume the 
whole process is a chain with k steps, </)j_i 2</)j, the last cf) will decay to two 
photons and cause 2k photons as the final products. In Appendix [XJ we give 
a general way to calculate the DOS for two-body decay chains with on-shell 
intermediate particles. In following subsections, we simply show some results 
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for different scenarios. 



2.1 DM + DM 27 and DM 27 

We first start with the simplest case where DM annihilates/decays directly 
to two photons. Since DM particles are almost stationary in the galactic 
frame, the photon DOS is a pure delta function, 

''-6{E,- — ) (2.1) 



N^dE^ V ^ 2 



Since = 2 the photon spectrum is simply twice the photon DOS. 



2.2 DM + DM 20 and DM 20 followed by ^ 27 

To compute the photon DOS when DM annihilates/decays to two bosons 
which subsequently decay to two photons each, we must first boost the DOS 
obtained in the previous subsection and convolute it with the appropriate 
DOS of DM annihilation/decay to two bosons as explained in Appendix 1X1 
Because the energy of the boson is always ^ and the direction of the 
photon in the boson rest frame is uniform, we can simply boost the delta 
function into the DM center of mass frame to obtain the photon DOS, 



1 dN^ _ 2 



(2.2) 



where the photon energy is between 




and is the boson mass. As expected the photon DOS is normalized to one 
and the photon spectrum is four times the photon DOS (A^-^ = 4). Notice 
that in the limit where m<^ = which case the two (p bosons are produced 

at rest in the DM center of mass frame, the photon DOS becomes a delta 
function and matches the photon DOS obtained in the previous subsection. 
This is easily understood since in this limit the bosons are stationary and 
decay to two photons each as in the previous subsection, only the delta 
function support and the DOS normalization change. 
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2.3 DM + DM 

and TT ^ 27 



20 and DM 20 followed by 



271 



When DM annihilates/decays to two bosons which decay to two vr bosons 
and finally decay to two photons, the photon DOS can again be computed by 
boosting the DOS obtained in the previous subsection and convoluting the 
result with the appropriate DOS of DM annihilation/decay to two bosons. As 
shown before, the (p DOS is a delta function and the n DOS is a step function. 
However the photon DOS gets complicated since, for a fixed photon energy, 
not all 7T with energy in the support of the step function can contribute. 
Indeed we have to calculate the minimum and maximum allowed vr energies 
which can generate the relevant photon energy and compute the photon DOS 
accordingly. Detailed computations are discussed in Appendix |Al Defining 
the boson masses as and the photon DOS is 



1 dN^ 



where 



Al 



In 
In 
In 



for ^^'^(E::^^^) <E^< ET'HET'' 



for Er'^f^^^''") < E^< E^^'^iEf'' 



for ET^'iEf'') < E^< ET^'iET'' 



(2.3) 



A 



1 - 



El 



El 




and = E^ — m^. The limits on the photon energy E.y are found using 



77' mm / 

^7 



Et, 




771 max 

^7 



T 



'1- 



El 



where El^'"" 



I Z?m: 



I Z?m: 



( Z?mi 



^min/max^^r 



is the minimum/maximum photon energy that can be generated by a vr boson 
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with minimum energy and similarly for £'™"^'™^^(^£'max^_ Again the photon 
DOS is canonically normalized while the photon spectrum is normalized such 
that the total number of photons is 8 (A'^ = 8) and it is shown in figure [T] 
for some given boson masses. 




Figure 1: Photon spectral distribution for DM + DM 20 and DM 20 
followed by ^ 27r and vr ^ 27 with M = 2000 GeV, = 400 GeV 
and = 0.14 GeV. The distributions peaks at 2E^/M = ^^/^^^^ = 
j^^max_|_ |p^ax|^^j^^gi/2^ whcrc 6 ^ 2.718; for the parameters here the peak is 
at 2E^/M ~ 0.581. 



There are three parts in DOS. Firstly, DOS increases with photon energy 
as a log function, then in the second part DOS is just a constant, does not de- 
pend on Ery, and in the third part DOS will decrease with E^ logarithmically. 
Since the parameters we take for the figure, the first part is out of region, 
the second part is corresponding to a straight line at low energy regime, and 
the joint point of second and third part will cause a kink at 2E^/M ^ 0.04 
since they are not smoothly connected to each other. 

As a consistency check, we can again match the DOS obtained here with 
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the DOS obtained in the previous subsections by taking different hmits on 
the mass ratios. Indeed, by talcing the hmits ^ = | or ^ = ^ (but not 
both) the DOS becomes a step function because, in such hmits, there are 
two particles produced at rest in the center of mass frame of the parent 
particle(s). Moreover in the limit where both ^ = ^ and -jj = \ the DOS 
becomes a delta function since, in that limit, DM annihilates/decays to two 
stationary (p bosons, and both decay to two stationary vr bosons, which 
finally decay to a total of eight photons, each of them taking | of the initial 
energy M. 

If we take some generic values for the two mass ratios — and the 
third part of the DOS, which decreases as — ln(£'^), dominates the DOS (it 
has the largest support). Also, adding more steps to the decay chain results 
in a softer photon DOS, a general trend which is intuitively expected since 
the total energy is distributed among a larger number of particles. 



2.4 DM + DM ^ iV0 and DM ^ iV0 followed by 



27 



Finally, we study DM annihilation/decay to N (p bosons which then decay 
to two photons each. The computation is done in the massless limit, i.e. 

= 0. This scenario occurs for example when two strongly-coupled bound 
states annihilate to very light pseudo Nambu-Goldstone bosons (pions in the 
QCD analogy) which then decay to two photons [3l]. In the massless limit, 
the boson DOS, assuming constant matrix element, can be computed by 
dimensional analysis and is given by 



dN^ 2{N -1){N ~2){M -2E^)^~^2E, 



MN-i 



(2.4) 



where Q < E^j, < Therefore the photon DOS is simply 



1 dN^ _ 2{N -1){M -2E^ 



\N~2 



dE^ 



(2.5) 



where < -E"^ < ^ (more detail is given in Appendix E])- The photon DOS 
satisfies = 1 and the photon spectrum is normalized to A'^ = 2A^. When 
N = 2 this result is equivalent to the scenario discussed in subsection 12.21 
equation fl2.2p . with = 0. The photon DOS for different A^ is shown in 
figure [21 One can easily see that the width grows with N. 
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Figure 2: Photon spectral distributions for DM + DM — ^ Ncj) and DM — ^ Ncj) 
followed by 27 in the limit M > m<^, for = 2, 3, 5, 10 and 20. All the 
distributions peak at 2NE^/M = E^/(E"^^^/A^) = 2. 



3 Photons from final states with charged par- 
ticles 

In this section, we analyze the irreducible photon background coming from 
charged particles, considering only the dominant process, i.e. FSR of photons 
or direct photon production from higher-order operators. 

For DM annihilation, FSR is the dominant process unless DM is a scalar 
boson or a Majorana fermion which annihilates directly to an electron-positron 
pair. For DM decay, FSR is the dominant process unless DM is a scalar boson 
which decays directly to an electron-positron pair. Indeed, in these specific 
cases FSR is small due to the large chiral suppression (in the s-wave ap- 
proximation for the annihilation scenarios), and might or might not be the 
dominant process according to the typical scale of the leading higher-order 
operators. 
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Here we will study two different modes which contribute to the photon 
DOS: DM annihilation/decay to one electron-positron pair and DM annihi- 
lation/decay to one boson pair which subsequently decay to one electron- 
positron pair each. 



3.1 DM + DM ^ e+ + + 7 and DM ^ e+ + + 7 



Final state radiation of photons 

For direct DM annihilation/decay to one electron-positron pair, the photon 
DOS from FSR of a single photon from the electron or the positron, in the 
collinear limit (a formula at leading order in the electron mass is given in 
Appendix [B]) , is simply given by 



and the DOS is shown in figure [3l Since there is soft photon Log divergence 
in the FSR, we choose to normalize respect to zeroth order approximation 
of a, i.e. o'DM+DM~^e++e- ioT annihilation and T£)M->e++e~ for decay. Notice 
here that the photon spectrum is the same than the photon DOS. 

Direct photon production from higher-order operators 

Photons from FSR may not always be the leading contribution to the spec- 
trum in some cases. The smallness of the electron mass might lead to large 
chiral suppression for FSR. So the leading process can be the one where 
photons are directly generated associated with e+e~. In low energy effective 
theory, one can write such process as higher order operators, assuming the 
particle intermediates the process to be heavy. We will show below that 
there is a low Mint regime where the higher-order operators dominate over 
FSR for the following scenarios: direct scalar boson or Majorana DM s-wave 
annihilation to one electron-positron pair or direct scalar boson DM decay 
to one electron-positron pair. 

For all these scenarios, named as scalar boson DM annihilation, Majorana 
DM annihilation and scalar boson DM decay, in the higher order operators, 
they share the common part as e^a^eF"'^ or eV'^eF"^. The only difference 
among them comes from the part of DM operators. But that part, under s- 
wave approximation, will only contribute different constant coefficients to the 
spectra of those scenarios. Since we are calculating DOS instead of spectrum, 
that difference will be exactly canceled by normalization. Thus, they share 
the same expression on DOS. 



1 dN.^ + (M - 2E^Y 



In 




(3.1) 
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Figure 3: Photon spectral distribution for DM + DM — e"^ + e +7 
and DM e+ + e" + 7 from FSR with M = 2000 GeV. The dis- 
tribution peaks at 2E^/M = E^/E^^^ = x where x is the solution of 
x{x'^ -2x + 2)/((l - x)(3x^ -Ax + 2)) = In [M'^{1 - x)/ml]; for the param- 
eters here the peak is at 2E^/M ^ 0.962. 



Within our assumptions, the photon DOS from the leading contact in- 
teractions given below in the limit of vanishing electron mass are all given 
by 



1 dN.^ _ 320 (M - 2E^)E^ 



dE^ M5 



(3.2) 



where E^^^ = M/2. The photon distribution for this spectrum is shown 
in figure HI The spectrum is independent of the specific case considered 
(scalar boson DM annihilation, Majorana DM annihilation or scalar boson 
DM decay). Here the limits on the photon energy are < E^ < ^ and DOS 
is normalized to 1 as well. Once again, the photon spectrum is the same than 
the photon DOS in this case. 
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Figure 4: Photon spectral distribution for DM + DM — + e~ + 7 and 
DM — > + e~ + 7 from the leading short range contact interaction for scalar 
boson or Majorana fermion S'-wave annihilation, or scalar boson decay, in the 
limit M > me. The distribution peaks at 2E^/M = = 5/6. 



Final state radiation versus direct photon production 

When scalar boson DM </> annihilates directly to one electron-positron pair, 
higher-order operators dominate for small Mint while FSR dominates for large 
Mint- Indeed, FSR is generated mainly by the dimension 6 operator 



while direct photon production is generated from dimension 8 operators 
through the following effective Lagrangian 



C,s D \-^d^{(j)^(j)){e^a''e)[aLF^, + bLp,,] + {L^R,e^e} (3.4) 



(3.3) 



int 
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where two-component spinor notation [oo] is used and the couphng constants 
{h,aL,bL,aR,bR} should naturally be order one numbers. A similar opera- 
tor has also been considered in Ref [31]. The relevant cross-sections are 



{(Tfsrv) ~ CFSRtt^^ln (^^j and (cTcfr) ~ Ccfra^^ where cfsr oc h"^ and 
Ceflf oc a| H — ■ include the appropriate order one coupling constants together 
with the 71 factors from the phase space integration. 

When Majorana DM x annihilates directly to one electron-positron pair, 
higher-order operators dominate for small Mint while FSR dominates for large 
Mint. Indeed, FSR is generated mainly by the dimension 6 operator 

£fsr D -^{x^(T''x){e^(r^e) + {L^R,e^e} (3.5) 

while direct photon production is generated from dimension 8 operators 
through the following effective Lagrangian 



\/AtTO! i J- 

jCcS D ^j^{x'cT''x)ie^(T''e)[aLF^, + fe^F^,] + {L^R,e^e} (3.6) 



int 



where two-component spinor notation [36] is used and the coupling constants 
{hL,aL,bL, hR,aR,bf{} should naturally be order one numbers. The rele- 
vant cross-sections are {(Tfsrv) ~ cfsr"]^ In (i^) and ((Tcfr) ~ Cefra^^^ 

where cfsr oc + and Cefr oc a\ + ■ ■ ■ include the appropriate order 
one coupling constants together with the tt factors from the phase space 
integration. 

Finally, when scalar boson DM decays directly to one electron-positron 
pair, higher-order operators dominate for small Mint while FSR dominates 
for large Mint. Indeed, FSR is generated mainly by the dimension 5 operator 

Cfsk D ^0(ee + ete-t) (3.7) 

while direct photon production is generated from dimension 7 operators 
through the following effective Lagrangian 

£eff 3 ^^9^0(eta^e) [a^F^. + bLF^,] + {L^R,e^e] (3.8) 

where two-component spinor notation [36] is used and the coupling constants 
{h,aL,bL,aR,bii} should naturally be order one numbers. This operator 
has also been considered in Ref [HI]. The relevant decay rates are Ffsr ~ 

cpsRa^TT?^ In ( ^ ) and Fefr ^ c^sa^ where cfsr oc h'^ and Cefr oc H 

int \ e / int 

12 



include the appropriate order one coupling constants together with the vr 
factors from the phase space integration. 

For all scenarios both processes (final state radiation and direct photon 
production) have comparable contributions when 

'^lO^GeV (3.9) 

where we assumed order one coupling constants and M 1 TeV. There- 
fore higher-order operators dominate over FSR for Mint ^ "while FSR 
dominates over higher-order operators for Mint ^ ^tit- Obviously for the 
scalar boson DM decay the decay rates are taken to be small enough such 
that DM is long-lived. This is possible if the operator coefficients are small, 
which occurs for example when the effective dimension of the scalar boson 
DM is higher. This is the case when the scalar boson DM is a composite 
field (like a glueball). In this case, equation 03.91) is still valid since the ef- 
fective dimensions of the operators for both final state radiation and direct 
photon production increase together. A more detailed analysis can be found 
in Appendix O 

3.2 DM + DM 20 and DM 20 followed by ^ 

e+ + + 7 

Here the computation of the photon DOS involves once more a convolution of 
the boosted (FSR or higher-order operators) photon DOS obtained in section 
13.11 with the appropriate DOS of DM annihilation/decay to two bosons as 
described in Appendix |X1 



M*, 
lilt 



M 



M 



In 



mi 



Final state radiation of photons 

Defining the boson mass as m^, the photon DOS from FSR of a single photon 
from the electrons or the positrons is thus 




(3.10) 
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where the function / is 




ln(x) + 2dilog(l 



x] 



M 



and the boson energy and momentum are — ^ 
respectively. The hmits on the photon energy are 



and 



M 
2 



(3.ii: 



4m 



2 

t 

Af2 



E, 



int 



Eris — 



nij 



and 



771 max 

^7 



E^ + 



The photon spectrum is simply twice the photon DOS [N^ = 2) which is 
shown in figure [5] for a given boson mass. 

Again it is possible to relate the FSR photon DOS of this subsection to 
the FSR photon DOS of the previous subsection by letting the boson mass 
approach half the parameter M, i.e. = M/2. This corresponds to DM 
annihilation/decay to two bosons at rest in the DM center of mass frame 
which subsequently decay to one electron-positron pair each with a single 
FSR photon. Notice that, in this limit, the {E'^ — m^)~2 prefactor in the 
photon DOS is very important for the matching to work. 



Direct photon production from higher-order operators 

As previously shown, photon production from scalar boson decay to electron- 
positron pair is not dominated by FSR when the typical scale of the interac- 
tions Mint is low. In that case, the photon DOS from higher-order operators 
is given by 

1 dN^ ^ I lor U < < — 2 — 

dE^ ~ I 5[(£^+|p^|)^-16(2i;^+2|p^|-3g^)g^] r E^-\p^\ . p . E^+\p^\ 

^ ^ { 3\p^\{E^+\p^\)i 2 ^ ^7 ^ 2 

(3.12) 



2 



where the boson energy and momentum are E,/, = ^ ci^nd — ^ v hP 
respectively. The photon DOS, which is shown in figure O is canonically 
normalized and the photon spectrum is simply twice the photon DOS (A'^ = 
2). Once more, it is possible to relate the photon DOS obtained in this 
subsection with the photon DOS obtained in the previous subsection for 
higher-order operators by taking the limit = y- 
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Figure 5: Photon spectral distribution for DM + DM 20 and DM 20 
followed by e+ + e" + 7 from FSR with M = 2000 GeV and = 400 
GeV. The distribution peaks at 2E^/M = E^/E"^^^ ~ 0.394. 



4 Photons from taus 

Tau leptons represent another interesting decay mode for dark matter anni- 
hilation/decay. Since baryon number is conserved in tau decays, and the tau 
mass is less than the sum of the proton and neutron masses, rrir < mp + rrin, 
tau decays include only (anti-) leptons and (anti-)mesons, with no (anti-) 
baryons. Dark matter that annihilates/decays preferentially to tau leptons 
is therefore not necessarily in conflict with stringent limits on the anti-proton 
flux in cosmic rays. However, tau leptons do provide an interesting source 
of photons since tau decays include a significant fraction of neutral pions, 
T —>■ X + 71^, that subsequently decay to photons, —>■ 27. 

The average number of tt^ produced in a single r decay is approximately 
(A'^o) ~ 0.51, giving an average of roughly one photon per r decay, N^/Nt ^ 
1. Most vr^ from r decays come from the hadronic one-prong decay modes 
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Figure 6: Photon spectral distribution for DM + DM 20 and DM 20 
followed by ^ e"*" + e~ +7 from the leading short range contact interaction 
with M = 2000 GeV and = 400 GeV, and neglecting the electron mass. 
The distribution peaks at 2E^/M = E^/E™^^ ~ 0.577. 



r~ — s> u^ + p^ z/T- + 7r~ + 7r'^ (branching fraction 25.4%) and Vr + (h 

Ut- + IT" + 271^ (branching fraction approximately 9% [SZ])- These two decay 
modes account for approximately 85% of all the neutral pions arising from 
r decay. The other main sources of neutral pions are the three-prong mode 
T~ —>■ Ur + 27r~ + TT"'" + 7r° (branching fraction 4.3%) and the one-prong 
mode T~ ^ Ut- + 7T~ + 3tt^ (branching fraction 1.1%), as well as continuum 
contributions. All branching fractions are taken from the PDG [3H]. Below 
we include only the dominant decays through the p and ai resonances. 

To obtain the photon DOS we first need to obtain the DOS of neutral 
pions. We have explicitly computed the contribution to the pion spectrum 
from the principal decay modes, with intermediate p and ai vector meson 
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resonances. The general r differential decay rate takes the form 



dV (X dHoiT ^ u-r + v) dml dUr 



mr] 



X |-M^,±(r + v)Vi;''{ml)M,{v ^ nir)]' (4.1) 

The 7r° DOS is obtained by removing a single ir^ from the final state phase 
space integration. Here V^'^{ml) is the vector meson propagator, Ai^,±{T —>■ 
Vr + v) denotes the matrix element for a r of helicity /i = ±| corresponding 
to right or left handed respectively, to decay to a vector meson, and M.y{v 
mr) denotes the matrix element for the vector meson to decay to n pions. 

First consider the decay through the p resonance. Obtaining the pion 
spectrum from the sequence of two-body cascades r— >-z/-|-p— s>z/-|-27r uses 
many of the techniques as used in the scalar cascades in Section [21 The 
major differences are first, nonconstant matrix elements, resulting from the 
nonzero spin of the intermediate p, and second, the large decay width of the 
p, which necessitates the use of a Breit-Wigner with a running width. We use 
a p mass and width of mo,p = 770 MeV and T^ p = 150 MeV. Next consider 
the decay through the ai resonance which differs from the p mode in that 
the final decay ai 3n is not two-body. The spectrum of the observed 
pion therefore requires additional integrations over the phase space of the 
unobserved pions in the final state. Following [371 EH], we use relatively 
simple parameterizations given in [IQ] for both the Oi Stt matrix element 
and running width. We use an ai mass and width of mo,a = 1.22 GeV and 
Fo.a = 420 MeV. 

The photon DOS may be obtained from the pion DOS by convolution, 
as before. We work in the coUinear limit M ^ m,- in which the components 
of the 7r°, and therefore photon, momentum transverse to the tau direction 
of motion in the origin annihilation/decay frame are irrelevant. The results 
for the normalized photon DOS under the assumption that the DM annihi- 
lates/decay to a single tau-anti-tau pair are well fit by the parameterized 
functional form 

= fiE,/E-ne-'^^^'/''^'^^ (4.2) 



where E"^^^ = M/2, and 



n=0 

3 

g{x) = J2cnx'' (4.3) 

n=l 
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For positive or right handed hehcity the best fit parameters are a+ = —0.192, 
= {3.90,-2.26,-15.59,32.96,-19.45}, c+,i = {8.94,-16.10,-19.45}, 
while for negative or left handed helicity the best fit parameters are a_ = 
-0.040, = {5.55,4.29,-10.94,0,0}, c_,i = {7.36,0,0}. The differences 
between the photon DOS arising from right- and left-handed taus are only 
very minor. It seems unlikely that a measurement of the photon spectrum 
alone could ever distinguish the helicity of tau's arising from DM annihi- 
lation/decay. The photon distribution arising from DM annihilation/decay 
directly to a single tau-anti-tau pair and averaging over tau helicity is shown 
in Fig. [71 In this case the photon spectrum is roughly twice the photon DOS, 



0.1 - 




2Ey/m 



Figure 7: Photon spectral distribution for DM + DM 2r and DM 2r 
followed by r — s> X + 7r° and vr*^ — > 27 in the limit M ^ m^. The distribution 
peaks at 2E^/M = E^/E;^^^ ^ 0.307. 

A'^ ~ 2, since there is roughly one photon per tau on average and two tau's 
per DM annihilation/decay. 

If we assume DM particles annihilate/decay to two intermediate bosons 
0, each of which subsequently decays to a tau-anti-tau pair, then we again 
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need to boost the DOS obtained above and convolve it with the appropriate 
DOS of DM annihilation/decay to two bosons, as explained in Appendix 1X1 
The final photon spectral distribution for this case in the collinear limits 
M ^ m0 ^ m,-, and averaging over tau helicities is shown in Fig. [HI The 
photon spectrum in this case is roughly four times the photon DOS, A'^^ ^ 4. 
We can see that the final photon DOS is softer compared to Fig. [7] due to 



0.1 r 




Figure 8: Photon spectral distribution for DM + DM 2(f) and DM 20 
followed by — 2r and t —>■ X + and 7r° 27 in the limit M ^ ^ 
rur. The distribution peaks at 2E^/M ~ 0.192. 

extra intermediate state. 

5 Photon spectra and flux 
5.1 Photon spectra 

Strickly speaking, everything we have computed up to now are photon DOS 
in different scenarios. The different photon DOS are superimposed in figure 
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[HI However, in the likely case where several scenarios add up, the differ- 




Figure 9: Photon spectral disbributions arising from different DM annihila- 
tion/decay scenarios for M = 2000 GeV, = 400 GeV and = 0.14 GeV. 
For DM annihilation/decay to intermediate (j) bosons the limit M ^ 
is presented with = 10. For DM annihilation/decay through intermediate 
r's the limits M ^ rrir or M ^ ^ m,- are presented. In each case the 
photon DOS is normalized to unit probability, except for FSR and Boosted 
FSR which are normalized with respect to the leading decay without FSR. 



ent spectra must be added with the appropriate weight. The full photon 
spectrum obtained is given by 

where the sum is over all allowed scenarios i, with photon spectrum 
and branching ratio Br(z). In comparing with the background of gamma ray 
photons, the only difference between DM annihilation and DM decay comes 
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from the power of the DM density profile in the appropriate equations as 
discussed in the next subsection. 



5.2 Photon flux 

An observation of gamma rays from dark matter annihilation/decay involves 
not only the spectrum but also the absolute magnitude. The photon flux 

- dA dt dn ^ ^ 

where dA and dQ are the detector area and solid angle elements, can be 



computed from an integral over the source along the line of sight (los) 
since gamma rays are not significantly attenuated on galactic length scales. 
For annihilation of a single species of DM particle that is its own anti-particle 
the spectral flux is given by 



(E,, n) = -t-^-^ \ pluir, n)dr (5.3) 



dE^ ^' ^^^Im J los 

where (av) is the phase space averaged annihilation cross section times ve- 
locity, pdm is the dark matter density, and dN^/dE^ is the photon spectrum 
with A'^^ = J dE^{dN^/ dE^) photons emitted per annihilation. If the DM 
is composed of distinct particle and anti-particle particle species with equal 
densities that can annihilate only through the particle-anti-particle channel, 
the flux fl5.3p should be multiplied by an additional factor of |. In terms of 
dimensional units the annihilation spectral flux (15.31) may be written 

cm-2 . s-1 ■ sr-1 VS X 10-26 cm3 ■ s-iy V ^dm / dE^ ^ ' 

(5.4) 

where 

( o.3 Gov. cm- J /o/^"'-"'''^ 

is a dimensionless order one factor that represents astrophysical parameters. 
For decay of a single species of DM the spectral flux is given by 

where F is the DM decay rate with A^ = J dE^{dN^/ dE^) photons emitted 
per decay. In terms of dimensional units the decay spectral flux may be 
written 

d^^{E^,n)/dE^ 9/ F \ /100GeV\ c/A^, 



cm ^ ■ s ^ ■ sr^ 



„/ F \ / 100 GeV\ dN.. , , 

^ ^ 6.26 X 10"^ — - -TT^M^) 5.7 
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where 

is a dimensionless order one factor characterizing the astrophysical parame- 
ters. With the spectral flux equations (15.31) and (15. 6p comparisons with the 
gamma ray background are straightforward. 



6 Conclusion 

High-energy photons are hardly deflected when they propagate through the 
galaxy. This is an advantage over charged particles, like electrons and positrons, 
which interact with the galactic magnetic field. Moreover, once charged par- 
ticles are accelerated, photons are always produced due to radiation. The 
different photon spectra one can obtained are good candidates to provide 
information on the process(es) generating the photons. In this paper, we 
studied the spectra of high-energy photons generated by various dark mat- 
ter annihilation/decay scenarios including: direct photon production from 
arbitrary two-body decay chains, final state radiation from charged parti- 
cles generated by dark matter annihilation/decay, direct photons production 
together with charged particles from higher-order operators, and the spe- 
cial case where photons are produced from taus generated by dark matter 
annihilation / decay. 

We noted also that, for processes which generate photons and leptons, 
effective field theory allows a comparison between the spectra from final 
state radiation and the ones from direct photon production due to higher- 
order operators. Interestingly, we found that, except for scalar dark matter 
annihilation or decay and Majorana dark matter annihilation, all the other 
dark matter scenarios are dominated by final state radiation. Moreover, 
for these exceptions (scalar boson dark matter annihilation and decay and 
Majorana dark matter annihilation), with a dark matter mass of (9(1 TeV) we 
found that direct photon production from higher-order operators dominates 
if the scale of the leading operator is lower than C(lOOOTeV). Otherwise 
final state radiation still dominates. Finally, it is also very interesting to 
see that the hardest spectrum among the spectra studied here (see figure [9]) 
comes from these exceptions, i.e. direct photon production from higher-order 
operators. 

Once the fiux of cosmic gamma rays is measured, an eventual dark matter 
signal could be compared with the different spectra presented here and a great 
deal of information on the nature of dark matter at the particle physics level 
could be deduced. 
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A Density of states 



In this appendix we review how the DOS is obtained from general consider- 
ations. The DOS for $ annihilation/decay to is simply given by 



1 dN 



1 d{av) 



N dE {av) dE 



and 



1 dN 

nIe 



1 dT 
fdE 



(A.l) 



respectively, and is thus normalized to one. For example, assuming constant 
averaged matrix element squared, the DOS in the $ center of mass frame for 
2^2 annihilation in the s-wave approximation and 1-^2 decay is simply 
obtained from the phase space. 



1 dN 

nIe 



d^p d^q 



EdE E 



X 



^5^^\P-'p-q) 



——6 (M-p-q) 



-\ -1 



(A.2) 



where P'* = (M, 0) with M = 2m$ for $ annihilation and M = m$ for $ 
decay. 

For $ decay in a boosted frame, the DOS follows from Lorentz covariance. 
Indeed N = J dE^ is a Lorentz scalar thus A^^^ = iv^°°***. Therefore one 
can rewrite the boosted DOS as 



1 dN'^°°'' 



N dE 



or 



CM 7 CM 



dE'-'^'dz 



1 dN^°°^^ 
N dE 



Ndz^^dE^^ 



dEdz 



CM 



d{E 



CM ^CM> 



dz 



CM 



d{E 



CM ^CM^ 



Ndz^^dE^^ 



Ndz^^dE^^' 



(A.3) 



Here z^^ = cos 6* = ]3$ ■ is the angle between the boosted $ and 
the unboosted (p. Since the DOS in the $ rest frame is uniform, then 
Ndz^dE<^^ = I NdE(^^ ■ Finally, the energy in the $ center of mass frame 
is related to the energy in the boosted frame by 



E 



E^ 



j^CM _^ ^CM 



_E'CM- 



rrij 



E'^ 



(A.4) 
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The bounds on z can be found from the bounds in the $ center of mass 
frame plus the physical constraint that —\<z<\. Boosting the DOS 
of the previous example in the frame where $ has four-momentum = 
E<^{\,f)<^^\ - m%/El), one gets 



1 dN^°°^^ 
N dE 



1 



4m2/m| - K 



where the energy is bounded by 

Ei^ 

E^ 

T 



4ml/ ml J 1 " ^l/^l 



Am1/m%\ll 



ml/ El 



(A.5) 



Convolution and matching of DOS 

For a two-body decay chain 0o ^ 20i to 4>k-i 20^ the DOS in the 0o rest 
frame can be found by iteration and is given by 



NdE. 



dE, 



fc-i ■ 



■■dEi 



NidEi N2dE2 ' ' ' NkdEk 



^jyBoost 



dEk^i- 



dNi 



k-l 



^jyBoost 



Nk-idEk-i NkdEk 
(A.6) 

where the bounds are complicated functions of the energies. For example, 
with the same assumptions as above, the two-body decay chain 00 — > 201 
followed by 0i — > 202 gives 

dN2 



N2dE2 



dEi 



dNi dNf°''^^ 
NidEi N2dE2 



[ dE,6 (e, - ^) 



- 4mi/m? /Bf 



m^ 



mo^yl — Amf/mly^l — Am\/m\ 
where the 02 energy is bounded by 



(A.7) 



771 max 
^2 



mo 
4 

mo 
4 



+ y'' 1 — 4mf /mg-y/ 1 — 4m2/mf 



Amj/mly 1 — 4m2/mf 



By taking limits where intermediate particles are created at rest in the 
center of mass frame of the parent particle, it is possible to match different 
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DOS. Indeed in the limit where rrii^i the decay chain is effectively 

cut by one step, with all the previous steps being unchanged, the i-th step 
being deleted and the subsequent steps being modified due to the energy 
redistribution. The DOS satisfies 



iim , ^ 



mi=2mi+i 

dNk^i 
Nk-idEk- 



(A.8) 



.,mi_i = 



-;mi=mi+i,...,mfc_i=mfc 



where the ith step is not included on the LHS. This allows us to check the 
different DOS formula. 



B FSR collinear divergence 

In the calculation of FSR, one would have to deal with the collinear diver- 
gence, which will show up when we take the massless limit for e'^e~ . Taking 
the electron mass into account, the three-body decay rate does not suffer 
from a collinear divergence and it can be easily computed using Dalitz co- 
ordinates [35]. Here we show the result for — > e"*" -|- e~ + 7, where FSR 
coming from the following vertex ^^0(ee + e^e^), one can get spectrum as 



dr ah^ml fml+{m^-2E. 



dE^ Sn^Ml, V "^0^7 



In 



ml 



2K-2E^^^^^, (B.l) 



E^ ^' 



Such spectrum has soft photon divergence, so we choose to normalize the 
spectrum respect to its zeroth order approximation on a, which corresponds 
to the process (f) + e~ , i.e. 



dN^ 1 dT 



'7 



where 



N^dE^ F0^e++e- dE,y 



(B.2) 



The first piece of the formula is from long distance contribution since this 
part is not sensitive to the details of the vertex, and it will be logarithmically 
diverge when we take m^ goes to zero. The second term is the finite piece 
when we take massless limit of e~^e~. This is the extra piece comparing 
to collinear approximation where we only keep long distance piece. In the 
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collinear limit (see equation fIC.ip ). after the last term of equation fIB.ip is 
dropped, one gets the photon DOS equation fl3.ip introduced in section [31 

Next, we consider the case where FSR gets one step of boosting, i.e. 
DM annihilate/decay to intermediate particle 0, then decays to e^e~ with 
photon from FSR. 



In a frame where has four-momentum p'^ = E^il.p^^l — m'l/E'^), the 

boosted DOS can be found from the general formula. However, since the 
photon is massless, the Jacobian simplifies greatly, leading to 



(B.3) 



The integral is easily done using the new variable w = - — cm"'/f2 =f) 

terms of which E^^ = wE^. The bounds on w can be found from the bounds 
on E^^ due to E^^{E^) and from the physical constraint — 1 < z < 1, giving 
two different regimes, 

. \ I Ami 

^ <w< mm < , ^ — ^ 1 - 



E, + ^El-ml {E, + ^El-ml^E,y m^ 

(B.4) 

as discussed in the text. Convoluting this boosted DOS with the two-body 
decay DOS one obtains the DOS for DM 20 followed by ^ e+ + e~ + 7 
mentioned above. 



C Higher-order operators 

In some specific scenarios, FSR suffers a large chiral suppression due to the 
small electron mass. In these cases, direct photon production from higher- 
order operators might dominate. In this Appendix we determine the neces- 
sary conditions for which direct photon production from higher-order oper- 
ators dominates over FSR. To reach the most general conclusions, effective 
field theory is used throughout the analysis. Therefore, the chirality rule, 
which states that chirality-violating operators must come with an overall 
mass term, is enforced. Moreover, to evaluate the appropriate FSR annihi- 
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lation cross-section/decay rate, the collinear limit is taken [12], 



dX. 



e++e~+7 

dx 



a 



TT 



[l-xf 



In 



X 



M2(l 



X) 



FSR 



vr 



In 



M2 



(C.l) 



where x = 2E^/M and X = {av) with M = 2mDM for DM annihilation or 
X = T with M = for DM decay respectively. Notice that for DM anni- 
hilation, DM is assumed to be almost at rest, thus the s-wave approximation 
can be taken. The lowest-dimension electron operators are given in table [H 
When the operator does not violate chirality, only the operators with e alone 
are considered since operators mixing e and e will lead to chirality-suppressed 
mixing terms. The photon can either appear in a covariant derivative, which 



0{x) 


dim 


(e+e-\0\ld) 




me(ee + e^e^) 


4 


me[{y-V+) + (xlx;)] 


4m^[(p_ -p+j - 77)2] 


ime{ee — e^e^) 


4 




4771^ ■ P+) + 777^] 




4 








4 








3 






a''(etCT''e) 


4 






i{D"e^ai'e - e*ai'D"e) 


4 








5 








5 






D^D^e^afe + eJai^D''D"e 


5 






iiD^D^e^ai'e - e^ai'D^D"e) 


5 







Table 1: Relevant electron-positron operators in the effective Lagrangian 
approach. The empty boxes correspond to operators which are not needed 
in the analysis. 



is already taken into account in the electron operators, or in the field strength 
tensor as in table [21 To simplify the analysis, the Lorentz indices are kept 



0{x) 


dim 


(7|O|0) 






2 




-(pX'ff'"' -PtK^"^"' -pX^''''' +p!;<ff'"'') 



Table 2: Relevant photon operators in the effective Lagrangian approach. 



free and are contracted with the appropriate tensors {g^y or e^j,Ap) only at 
the end of the analysis. In this way, one does not have to deal with the com- 
plete set of operators at this stage of the analysis (for example, the operator 
corresponding to the dual field strength tensor can be forgotten). 
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Scalar boson DM annihilation to e+ + e +7 

The relevant lowest- dimension scalar boson operators are given in table [3l 
Combining operators of table [3] with operators of table [H the leading op- 



0{x) 


dim 


{o\o\<p^) 




(\Mn 

\ / s— wave 




2 


1 


1 


1 




3 


-i{p': + p^) 


{pt+P'^W* +P''') 






3 


{p'i - p^) 


{p^-P^'M -P"') 






Table 3: Relevant scalar boson DM operators for DM annihilation. 



erator for FSR seems to be constructed from 9^(0^0) and (e^a^e). It is 
a dimension 6 operator and it is not explicitly chirality-suppressed. How- 
ever, in the s-wave approximation, p'^ = {m^^'Pe) and p^ = (m^, —pe), thus 
|(e+e~|(e^o"°e)|0)P = 2mg which is chirality-suppressed. Therefore, the lead- 
ing FSR operators are dimension 6 but are chirality-suppressed. For example, 
one such operator can be constructed out of 0V ^^^d me(ee + e^e'^). The re- 
lated FSR cross-section is 



2 



Am J 



ml _| 



(C.2) 



where Cfsr includes the operator coupling constant and the vr factors from 
the phase space integration. 

For direct photon production from higher-order operators to dominate, 
the operators should not be chirality-suppressed and/or should be of lower 
dimension. Combining again operators of table [3] with operators of tables 
[U and [2] such that one photon can be created in the final state, the leading 
operators for direct photon production are dimension 8 and are not chirality- 
suppressed. The operator made out of d^{(l)^(j)), {e^aye) and F^y is an ex- 
ample. Following the general rules of effective field theory, a factor of a 
should be put in front of this type of operators. The related direct photon 
production cross-section is 

{2m^f 

((Jcfft') ~ CcffQ g (C.3) 

where Cefr includes the operator coupling constant and the tt factors from the 
phase space integratioiu. 



"'^cfsr and Ces have the same phase space factors since both are obtained from a three- 
particle final state. 
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It is now straightforward to compare FSR and direct photon production 
cross-sections, 



cfsr mlM-L 



rat 



(C.4) 



which leads to equation fl3.9l) of section [HI Direct photon production from 
higher-order operators will therefore dominate over FSR when Mint ^ 
(see equation fl3.9p ). In the remaining subsections, the analysis is more con- 
cise since it closely follows what has been done here. 



Majorana DM annihilation to e+ + e +7 

The relevant lowest-dimension Majorana operators are given in tablelH Combining 



0(.t) 


dim 


(0|O|xx) 




\ / s— wave 




4 


2m^[{xiX2) + [y\yl)] 


■l'n>[(Pi ■P2) - ml] 







4 


2im^[{xiX2) - (y\yl)] 


'^ml[{pl ■p2) + ml] 




m^{x'^''-'X + X^^'"'X^) 


4 











im^{x'y'"'X-X^^>"'X^) 


4 











X^iri^X 


3 


[ylafx^) - [yla'^xi] 


pi'p'/ + p^j/;' - ■ P2) + mlgfi^' 


2mlSl'°S^^'" 




4 


—i{Pi + P2) ^ i^i^ 


{Pi+P2){Pi +P2) X ibid 


8m^5''''5'''"r"(5"'" 


i{afx^a>'X-X^^''d"x) 


4 


{P1- P2) X {- ^ +} 


[Pi - P2) W - P2) X {+™' ^ -™'} 






Table 4: Relevant Majorana DM operators for DM annihilation. Notice that 



since x is a Majorana fermion. 



operators of table H] with operators of table [T] the leading order operators for 
FSR are dimension 6 operators of the form {x^(T'^x){^^'^fi^)- Again, these op- 
erators are chirality-suppressed in the s-wave approximation, thus the related 
FSR cross-section is 



2 



4m 



(C.5) 



where cfsr includes the operator coupling constant and the vr factors from 
the phase space integration. 

Combining operators of table H] with operators of tables [T] and [21 the 
leading operators for direct photon production are dimension 8 and are not 
chirality-suppressed. The operator made out of {x^(T^x)y (e^cx^'e) and F^i, is 
an example. The related direct photon production cross-section is 

(cTeflfW) ~ Ceffa^^^ (C.6) 
int 
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where Ceff includes the operator couphng constant and the vr factors from the 
phase space integration. 

Comparing FSR and direct photon production cross-sections leads to the 
equivalent of equation ( 1C.4I) with the same overall conclusions as scalar boson 
DM annihilation. 



Dirac DM annihilation to e+ + e +7 

The relevant lowest-dimension Dirac operators are given in table Combining 



0{x) 


dim 


(0|O|rm) 


\MV- 


(imp) 

\ / s— wave 




3 










3 


{ylai'x) + (■!/tCT''x,) 


P^Pt' + Kp^' - 9''^' {P ■ P.) - m^S^^' 





Table 5: Relevant Majorana DM operators for DM annihilation. 



operators of table [5] with operators of table [H the leading order operators 
for FSR are dimension 6 operators of the form [(77 "'"o"'^?]) — (f/"'"(T'^f/)](e^a^e). 
However, these operators are not chirality-suppressed in the s-wave approxi- 
mation, thus the related FSR cross-section is 



{apsRv) ^ CFSRavTT In 

^'^int 



4m: 



(C.7) 



where Cfsr includes the operator coupling constant and the vr factors from 
the phase space integration. 

Combining operators of table with operators of tables [T] and [21 the 
leading operators for direct photon production are dimension 8 and are not 
chirality-suppressed, as in the Majorana DM annihilation case. The operator 
made out of [(?7^o"'^?7) -|- {fj'^a'^f])], {e^a'^e) and F^i, is an example. The related 
direct photon production cross-section is 

(defrf) ~ Ceflpa g (C.8) 

int 

where Ceff includes the operator coupling constant and the vr factors from the 
phase space integration. 

Comparing FSR and direct photon production cross-sections leads to 



CFSR Mil r4m2 



H 

ml 



> ((Tefft^) (C.9) 



(aFSR.)-(aeff.)^^ln 

and, for order one coefficients, FSR always dominates over direct photon 
production from higher-order operators. 
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Gauge boson DM annihilation to e+ + e +7 

The relevant lowest- dimension gauge boson operators are given in table El 
Combining operators of table [6] with operators of table [H the leading order 



0{x) 


dim 


{Q\0\XX) 




{\M\A 

\ ' s— wave 




4 


-[Ker-pr65')(P^4-f^4) + {l-2}] 








5 




[pl + pi)(Pi +pf)-x ibid 


4m2^-(5«"i»'" X ibid 



Table 6: Relevant gauge boson DM operators for DM annihilation. 



e- 



operators for FSR are dimension 8 operators of the form X^^Xy\i{D^ e^o ^ 
e^cjy^D^e). However, these operators are not chirality-suppressed in the s- 
wave approximation, thus the related FSR cross-section is 



((Tfsr^^) ~ cpsRa— TTg — In 

int 



ml 



(C.IO) 



where cfsr includes the operator coupling constant and the tt factors from 
the phase space integration. 

Combining operators of table E] with operators of tables [T] and [21 the 
leading operators for direct photon production are dimension 8 and are not 
chirality-suppressed. The operator made out of X^^Xy\ and iiD^e^Gp^e — 
e^a^D^e) is an example. The related direct photon production cross-section 
is 

(aefff) ^ Ccfrtt— -g — (C.ll) 

int 

where Cefr includes the operator coupling constant and the vr factors from the 
phase space integration. 

Comparing FSR and direct photon production cross-sections leads to 



{crFSRv) ~ (a-cfT^^)^^ln 

Ceff 



ml 



> 



{a,sv) (C.12) 



and, for order one coefficients, FSR always dominates over direct photon 
production from higher-order operators, although only slightly. 



Scalar boson DM decay to e+ + e +7 

This case follows closely the scalar boson DM annihilation case. The relevant 
lowest-dimension scalar boson operators are given in table [71 Combining 
operators of table [7] with operators of table [H the leading order operators for 
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0{x) 


dim 


(O|O|0) 






s— wave 




1 


1 


1 


1 




2 









Table 7: Relevant scalar boson DM operators for DM decay. 



FSR are dimension 5 operators of the form me(f>{ee + e^e^). These operators 
are chirality-suppressed and thus the related FSR decay rate is 



FSR 



2 H 



(C.13) 



where CpsR includes the operator coupling constant and the vr factors from 
the phase space integration. 

Combining operators of table [7] with operators of tables [1] and O the 
leading operators for direct photon production are dimension 7 and are not 
chirality-suppressed. The operator made out of d^(f>, {e^a^e) and F^y is an 



example. The related direct photon production decay rate is 

,7 



efr 



Cefra 



ML 



(C.14) 



where c^s includes the operator coupling constant and the vr factors from the 
phase space integration. 

Comparing FSR and direct photon production decay rates leads to the 
equivalent of equation (lC.4p for decay rates, 



FSR 



r 



cff 



cfsr mlML 



4m ^ 



with the same overall conclusions as scalar boson DM annihilation. 



(C.15) 



Abelian gauge boson DM decay to e+ + e +7 

The relevant lowest- dimension gauge boson operators are given in table El 
Combining operators of table [HI with operators of table [Tj the leading order 
operators for FSR are dimension 6 operators of the form X ^jj{e^ a'^ e) . How- 
ever, these operators are not chirality-suppressed and thus the related FSR 
decay rate is 



TfSR ~ CFSRtt 



M 



X 



In 



int 



m 



X 



mt 



(C.16) 
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Oix) 


dim 


(0|O|X) 


\M\^ 


(\Mn 

\ / s— wave 




2 










3 


—ip'^ X ibid 


p^p^ X ibid 


^I^AO^A'O X ibid 



Table 8: Relevant gauge boson DM operators for DM decay. 



where CpsR includes the operator coupling constant and the tt factors from 
the phase space integration. 

Combining operators of table M with operators of tables [U and Ej the 
leading operators for direct photon production are dimension 6 and are 
not chirality-suppressed. The operator made out of X^'^ and i{Di,e^a^e — 
e'^afj^D^e) is an example. The related direct photon production decay rate is 



5 



eff^Cefra^^ (C.17) 



where Cefr includes the operator coupling constant and the vr factors from the 
phase space integration. 

Comparing FSR and direct photon production cross-sections leads to 



CpSR , 

FSR ~ f cff In 

CefT 



> Teff (C.18) 



and, for order one coefficients, FSR always dominates over direct photon 
production from higher-order operators, although only slightly. 



DOS for direct photon production from higher-order 
operators: An example for Majorana DM annihilation 

When direct photon production from higher-order operators dominates over 
FSR, the photon DOS can be obtained from the most general effective La- 
grangian. Here we study the case of Majorana DM annihilation to electron- 
positron pair, X + X — > e"*" + e~ + 7. The other cases where direct photon 
production from higher-order operators dominates over FSR are basically 
equivalent. 

The lowest-dimension operators relevant to x + X ~^ + e~ + 'y are 
effective dimension 8 operators due to the chirality rule, thus only effective 
mass dimension 8 operators will be considered. The highest possible effective 
mass dimension for the DM (x) operators is therefore 4. The photon can 
either appear in a covariant derivative or in the field strength tensor. 
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With the help of tables [H [2] and HI the minimal set (in the sense that 
any operator relevant for x + X~*'6^ + e~ + 7 can be rewritten as a linear 
combination of the operators in the minimal set) of operators of the mass 
dimension 8 effective Lagrangian for the process X + X'~*6"'~ + e~+7 can be 
found. 

From the electron operators without covariant derivatives, only the one 
with mass dimension 3 is relevant, leading to {x^a'^x)i^^^'^^)F'^''- From the 
electron operators with one covariant derivative, only i{x^d'^x)9''{D^e^d''^e — 
e^a'^D^e) survives (there are two different ways of building this operator) 
since all other operators are either zero to lowest order in me from the equa- 
tions of motion (the a matrix and the covariant derivative are forced to be 
contracted together) or vanish in the s-wave. Finally there are only two oper- 
ators that can be built from the electron operators with two covariant deriva- 
tives, which are {x^a^'x){DPD^e^a''e+e^a'' DPD^e) axidi{x^a>'x){DPD^e^a''e- 
e^a^DPD^e). 

All operators have four Lorentz indices and must therefore be contracted 
with g^^^gxp, g/^xgup, g^^pgux and e^^Ap- Using the equations of motion, this 
fact leads to an even smaller minimal set (to lowest order in mg), since 
g^yl = a^ay + 2ia^y and [D^, D^] = — iv47raF^j,. Indeed one has 

gxpie^^^DPD^e) = e^a^iapaxD^D^ + g(JpxFP^)e ~ F^^'-term + 0(me) 
g^p^e^a^'DPD^e) = g^pC^ a" {D^DP + [DP, D^])e F^p -term + 0{m,) 
g^xie^a^DPO^e) ~ 0{m,) 

ep^Xpie^cT^OPD^e) = p,xp{e^ a'' [Dp , D^]e) F^p -term 

thus eliminating all operators with two covariant derivatives. Moreover, the 
operator with one covariant derivative can be rewritten in terms of two co- 
variant derivatives as i{x^ (jPx){F)'^ (^'^ D^e — D^e^a'^ DPe), leading to 



gxp{DPe^(7''D^e - D^e^cx'DPe) = 

g^pi^DPe^a^D^e - D^e^a^DPe) ~ 

g^x{D''e^a''D^e- D^e^a'^DPe) ~ 

epuXp{DPe^rD^e - D^e^a^DPe) = 





2ep,Xp{DPe^<y"D\) 
2ep,Xp[dP{e^a-D^e) - {e^a^'DPD^)] 
ep^XpO^ie^a^D^e) + F^'^-term 



therefore eliminating three of the four possible operators. Finally, for the 
operators with F^^ , only two operators survive since the field strength tensor 
is antisymmetric. 
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Then the minimal set consists of {x^a^x){^'''^'^^)^tJ.u, (x^'^^x)(cV''e)F^,^ 
and iei,^xpd''{x^a''x){e^^''D^e). In the s-wave approximation, the last term 
also vanishes, thus leading to only two operators relevant for the mass di- 
mension 8 effective Lagrangian of x + X c"^ + e_ + 7 (at lowest order in 



^cff = ^^{x^^^xKe^cT''e)[aLF,, + BlF,,] + {L^R,e^e} (C.19) 



ML 



where the coupling constants are assumed to be order one numbers. Since 
the operators do not interfere (they couple the electrons to different photon 
states), the probability is simply given by 



\M\' = -32Mal + bl + al + bl)mlM.J 

X [{Pj ■ P+){Pi ■ P-) - E^E+{p., ■ p_) - E.,E_{p^ ■ 
= 6i7ia{al + bl + al + bl)ml{m^-E,)Mr^ 

X {2ml ~ '^rn-xE^ - 4m.^E+ + E^ + 2E^E+ + 2El) 



where p_ = -p^ - p+, E^ = J E?^ + El + 2zE^E+ = 2m^ - E^ - E^ 



z = ^"^x 2m;^E-y^ 2mxE++E~,E+ — Ey < E^ < m^. Here z = COS 9 and 9 

is the angle between the positron and the photon. In the vanishing electron 
mass limit the annihilation cross-section is thus 

(C.20) 

with < Ey < m^, which gives the DOS mentioned in the text. Notice that 
the annihilation cross-section vanishes like the cube of the photon energy, 
E^, as Ey 0. 
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